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Abstract 

We perform a systematic analysis of the density dependence of the nuclear symmetry energy 
within the microscopic Brueckner-Hartree-Fock (BHF) approach using the realistic Argonne V18 
nucleon-nucleon potential plus a phenomenological three body force of Urbana type. Our results 
are compared thoroughly to those arising from several Skyrme and relativistic effective models. 
The values of the parameters characterizing the BHF equation of state of isospin asymmetric 
nuclear matter fall within the trends predicted by those models and are compatible with recent 
constraints coming from heavy ion collisions, giant monopole resonances or isobaric analog states. 
In particular we find a value of the slope parameter L = 66.5 MeV, compatible with recent exper- 
imental constraints from isospin diffusion, L = 88 ± 25 MeV. The correlation between the neutron 
skin thickness of neutron-rich isotopes and the slope, L, and curvature, K sym , parameters of the 
symmetry energy is studied. Our BHF results are in very good agreement with the correlations 
already predicted by other authors using non-relativistic and relativistic effective models. The 
correlations of these two parameters and the neutron skin thickness with the transition density 
from non-uniform to /3-stable matter in neutron stars are also analyzed. Our results confirm that 
there is an inverse correlation between the neutron skin thickness and the transition density. 

PACS numbers: 21.65.Cd; 21.65.Ef; 21.65. Mn 
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I. INTRODUCTION 



A well-grounded understanding of the properties of isospin-rich nuclear matter is a nec- 
essary ingredient for the advancement of both nuclear physics and astrophysics. Isospin 
asymmetric nuclear matter is present in nuclei, especially in those far away from the sta- 
bility line, and in astrophysical systems, particularly in neutron stars. A major scientific 
effort is being carried out at an international level to study experimentally the properties of 
asymmetric nuclear systems. Laboratory measurements, such as those running or planned to 
run in the existing or the next-generation, radioactive ion beam facilities at CSR (China), 
FAIR (Germany), RIKEN (Japan), SPIRAL2 / G ANIL (France) and the upcoming FRIB 
USA), can probe the behavior of the symmetry energy close and above saturation density 
1|. Moreover, the 208 Pb Radius Experiment (PREX), scheduled to run at JLab in early 
2010, should provide a very accurate measurement of the neutron skin thickness in lead 

n 

via parity violating electron scattering [2|. Astrophysical observations of compact objects 
are also a window into both the bulk and the microscopic properties of nuclear matter at 
extreme isospin asymmetries sj]. The symmetry energy determines to a large extent the 
composition of /5-stable matter and therefore the structure and mass of a neutron star {4J]. 

The empirical knowledge gathered from all these sources should be helpful in identifying 
the major issues arising when the isospin content of nuclear systems is altered. Reliable 
theoretical investigations of neutron-rich (and possibly proton-rich) systems are therefore 
called for. Phenomenological approaches, either relativistic or non-relativistic, are based on 
effective interactions that are frequently built in order to reproduce the properties of nuclei 
{5]. Since many of such interactions are built to describe systems close to the symmetric 
case, predictions at high asymmetries should be taken with care. A priori, the starting point 
of microscopic approaches appears to be safer: realistic nucleon-nucleon (NN) interactions 
reproduce the scattering and bound state properties of the free two-nucleon system and 
include naturally an isospin dependence jf| . The in-medium correlations are then built using 
many-body techniques that microscopically account for isospin asymmetry effects such as, for 
instance, the difference in the Pauli blocking factors of neutrons and protons in asymmetric 
systems 

In practical applications, phenomenological approaches can be significantly improved in 
the isospin asymmetric case by using, as input, microscopically based predictions. The 
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Skyrme interactions of the Lyon group j^l, 0, 10] for instance, reproduce a neutron mat 
ter equation of state (EoS) based on microscopic variational calculations Q and are, 
therefore, able to predict reasonable properties for compact stars 131 ]. 

Even when these restrictions are taken into account, however, some of the properties of 
asymmetric nuclear matter appear to be rather unconstrained. In particular, the density 
dependence of the symmetry energy is still an important source of uncertainties. Different 
approaches predict similar asymmetry properties close to saturation, but strongly diverge 
for densities either above or below the saturation point. We shall try to give a quantitative 
prediction for the density dependences arising from a microscopic perspective. Let us note 
that the situation in symmetric nuclear matter is quite different: the saturation density, 
binding energy and incompressibility are relatively well settled from an empirical point of 
view. However, in microscopic calculations [l4, [l6], the prediction of these saturation 
properties is strongly influenced by three-body forces (TBF), mainly concerning the deter- 
mination of the saturation density that can be easily off by 40% in the absence of TBF. 
Other microscopic calculations using also realistic interactions have been recently reported 
in the framework of Dirac-Brueckner-Hartree-Fock [nj, [3] 

Values for the properties of asymmetric nuclear matter can be obtained from various 
analyses of experimental data^ncluding isospin diffusion measurements [l| , g iant resonances 
[3], isobaric analog states {20I or meson production (pions [211 ]. kaons jzj) in heavy ion 
collisions. Another important tool to determine empirically these properties are the existing 
correlations between different quantities in bulk matter and finite nuclei. The Typel-Brown 
correlation, for instance, is a linear relation between the density derivative of the neutron 
matter EoS at 0.1 fm -3 and the neutron skin thickness of 208 Pb that has been theoretically 
tested with different Skyrme parameter sets j^J and relativistic Hartree models j^j]. Accu- 
rate measurements of neutron skin thicknesses, via future parity violating experiments 



or by means of existing antiprotonic atoms data 25 



26| . are thus helpful in determining 



the bulk properties of nuclear systems. Other linear correlations, such as those relating the 
208 Pb skin thickness and the liquid-to-solid transition density in neutron stars 221] , or power 
aw correlations, such as the relation between the radius of a neutron star mass and the EoS 



281 ] , have also been observed. There is so much dispersion on the results of these correlations 



obtained with phenomenological approaches, that fully microscopic calculations, as the one 
performed in this paper, are needed. 
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In the present work, we compute the density dependence of the symmetry energy and 



physical quantities directly related wit 
interaction, namely the Argonne VI 8 



i its slope and curvature obtained from a realistic 
291 ] plus a TBF of the Urbana type, in the frame- 
work of the Brueckner-Hartree-Fock (BHF) approximation. After a brief discussion of the 
parametrization of the energy density and symmetry energy, and a short description of the 
BHF approach, we discuss and compare extensively our results with those obtained with 
several Skyrme forces and relativistic effective models. We pay particular attention to the 
trends established by these calculations and analyze different correlations arising from them. 
We conclude the discussion by summarizing the more important results. 

II. ASYMMETRIC NUCLEAR MATTER 

Assuming charge symmetry for nuclear forces, the energy per particle of asymmetric 
nuclear matter can be expanded on the isospin asymmetry parameter, (3 = (N — Z)/ (N + 
Z) = (p n — Pp)/p, around the values of symmetric (/3 = 0) nuclear matter, in terms of even 
powers of j3 as 

0) = E SNM (p) + S 2 (p)(3 2 + S,(p)(3 4 + 0(6) , (1) 
where Esnm(p) is the energy per particle of symmetric matter, S 2 (p) is identified (excluding 



2d 



surface contributions 
formula 

ld 2 E/A 

o 2 \P) = 

and 



261 ]) with the usual symmetry energy in the semi-empirical mass 



aw = - 2 -^r „_„ , (2) 



(3=0 



1 d A E/A 

It is well known, however, that the dominant dependence of the energy per particle of 
asymmetric nuclear matter on (3, is essentially quadratic [3,0, 31 1 . Therefore, contributions 



from 5*4 and higher other terms can be neglected, and one can, in good approximation, 
estimate the symmetry energy simply from the two extreme cases of both pure neutron 
matter and symmetric nuclear matter according to 

S 2 (p)~|(p,l)-f(p,0). (4) 

In Fig. [1] we show the density dependences of the coefficients S 2 and S4 obtained in our BHF 
calculation (left panel), together with the results predicted by the Skyrme force SLy230a 



middle panel) and the relativistic mean field model NL3 (right panel). For the three models, 



32j l in the whole range of densities explored 



5*2 is an increasing function of the density 
(0 — 0.3 fm -3 ). The rate at which S2 increases is, however, substantially different for each 
of them: while NL3 predicts a steep, almost linear increase, SLy230a shows a substantial 
down bending above saturation. The BHF results appear to be somewhere in the middle of 
the two behaviors. Note that in the three expected, the coefficient S4 is very small 

(below 0.5 MeV in the BHF case, and below 1 — 2 MeV in the case of the Skyrme force and 
the NL3 model) in the density region considered. 

It is common to characterize the density dependence of the energy per particle of symmet- 
ric matter around the saturation density p in terms of a few bulk parameters by expanding 
it in a Taylor series around p , 

«)^0 + f(^) 2 + f(^) 3 + O (4). (5, 

The coefficients denote, respectively, the energy per particle, the incompressibility coefficient 
and the third derivative of symmetric matter at saturation, 

E = E SNM (p = po) , K = 9p — , Q = 27pg — . (6) 

Op 1 P=Po Op 6 P=Po 

Similarly, the behaviour of the symmetry energy around saturation can be also charac- 
terized in terms of a few bulk parameters, 

where E sym is the value of the symmetry energy at saturation and the quantities L, K sym and 
Q S ym are related to its slope, curvature and third derivative, respectively, at such density, 

L ~ 3p °^p~ P=P0 ' ^ ~ 9P °^~ P=P0 ' Qsym ~ 27P °^~ P=P0 ■ (8) 

Combining the expansions (PQ), (J5J) and (J7j), one can predict the existence of a saturation 
density, i.e., a zero pressure condition, for a given asymmetry, and rewrite the energy per 
particle of asymmetric matter around the new saturation density po(P) ~ p (l — 3(L / K ) (3 2 ) 
in a form similar to Eq. (JSJ), 
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where the coefficients Eo({3), K (j3) and Qo{(3) that characterize the energy per particle, the 
incompressibility coefficient and the third derivative at po((3) for a given asymmetry (3 read 
up to second order 

E Q (/3) = E + E sym (3 2 + O(A) 

K ((3) =K + [K sym -6L- Q±L) (3 2 + (9(4) (10) 
QM = Qo + (Qsym - 9L^j (3 2 + 0(4) . 

Fig. [2] shows the saturation density (left panel), energy per particle (middle panel) and 
incompressibility (right panel) as a function of /3 2 , up to a value of (3 ~ 0.6. For {3 = one 
recovers the results of symmetric nuclear matter, then as (3 increases the saturation density, 
the binding energy and the incompressibility decrease. These behaviors are rather intuitive 
and a direct consequence of Eq. fTTOT) and the specific values of E sym , L, K sym , K and Q 
at the saturation density of symmetric nuclear matter. Note that the BHF results are well 
reproduced by the quadratic expansion on (3 in the range of asymmetries considered. 

In the following, before presenting our results, we shortly review the BHF approach of 
asymmetric nuclear matter and provide a few details on the Skyrme forces and the relativistic 
models considered. 



A. The BHF approach of asymmetric matter 

The BHF approach of asymmetric nuclear matter [?], Q] starts with the construction of 
all the G matrices describing the effective interaction between two nucleons in the presence 
of a surrounding medium. They are obtained by solving the well known Bethe-Goldstone 
equation 

GtiT2;t 3 T4 (k>) — ^ / tiT2;t 3 T4 + / , V T1TrnT . ■ ■ —G TiTj - T3T4 (uj) (H) 

7f ' ^ - ( ; - ( j ■ III 

where r = n,p indicates the isospin projection of the two nucleons in the initial, intermediate 
and final states, V denotes the bare NN interaction, Q TiTj the Pauli operator that allows 
only intermediate states compatible with the Pauli principle, and u, the so-called starting 
energy, corresponds to the sum of non-relativistic energies of the interacting nucleons. The 
single-particle energy e T of a nucleon with momentum k is given by 

e T (k) = ^ + Re[U T (k)} , (12) 



where the single-particle potential U T (k) represents the mean field "felt" by a nucleon due 
to its interaction with the other nucleons of the medium. In the BHF approximation, U(k) 
is calculated through the "on-shell energy" G-matrix, and is given by 

U T (k) = J2 E 0' I G TT , TT r(u = e T (k) + e T ,{k'))\kk') A (13) 

where the sum runs over all neutron and proton occupied states and where the matrix ele- 
ments are properly antisymmetrized. We note here that the so-called continuous prescription 
has been adopted for the single-particle potential when solving the Bethe-Goldstone equa- 



tion. As shown in Refs. 



34 



35j , the contribution to the energy per particle from three- hole 
line diagrams is minimized in this prescription. Once a self-consistent solution of Eqs. (ITT|) 
and (fTBl is achieved, the energy per particle can be calculated as 

!<*/»= E (i4) 



The BHF calculation carried out in this work uses the realistic Argonne V18 (Avl8) 29] 
nucleon-nucleon interaction supplemented with a three-body force of Urbana type which for 
the use in BHF calculations was reduced to a two-body density dependent force by averaging 
over the third nucleon in the medium {36]. This three-body force contains two parameters 
that are fixed by requiring that the BHF calculation reproduces the energy and saturation 
density of symmetric nuclear matter (see Refs. 14j, |l5j, [16( for a recent analysis of the use of 
three-body forces in nuclear and neutron matter). Note that the Avl8 interaction contains 
terms that break explicitly isospin symmetry. Therefore, in principle, we should consider 
also odd powers of /3 in the expansion (pQ) in our Brueckner calculation. However, we have 
neglected such terms since, as shown by Muther et al. in Ref. [6j], the effects of isospin 
symmetry breaking in the symmetry energy are almost negligible (less than 0.5 MeV for a 
wide range of NN interactions). 



B. Phenomenological Models 

Effective nucleon-nucleon interactions of Skyrme type are very popular in nuclear struc- 



ture calculations (see Ref. [5] for a recent review). The zero-range nature of this phenomeno- 
logical NN interaction allows for a very efficient implementation of mean field calculations, 
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in both finite nuclei and extended nuclear matter where one can get simple analytical ex- 
pressions. In fact, the main advantage of these forces comes from their analytical character, 
which make them very useful to get a physical insight into problems where the fully mi- 
croscopic calculations are very difficult to implement. By construction, most of the Skyrme 
forces used in the literature are well behaved around the saturation density of nuclear matter 
and for moderate isospin asymmetries. However, not all Skyrme parameters are completely 
well determined through the fits of given sets of data and only certain combinations of the 
parameters are really empirically determined. This leads to a scenario where, for instance, 
different Skyrme forces produce similar equations of state for symmetric nuclear matter 
but very different results for neutron matter. Recently, an extensive and systematic study 
has tested the capabilities of almost 90 existing Skyrme parametrizations to provide good 
neutron-star properties. It was found that only twenty seven of these forces passed the 
restrictive tests imposed, the key property being the behavior of the symmetry energy with 
density fl3l ]. In this work, we have only considered the forces that passed the tests imposed 
by Stone et al. , in Ref . 13| . Particular numerical results concerning some forces of the Lyon 
group, the Ski family (adjusted to isotopic shifts in the lead region) [37|, |38( and the early 
SIII and SV |39j parametrizations will be discussed in Table [I] 

Regarding the relativistic effective approaches, we have considered two different types of 
models: (i) non-linear Walecka models (NLWM) with constant couplings, and (ii) density 
dependent hadronic models (DDH) with density dependent coupling constants. Within the 
first type, in particular, we have considered the models NL3 [40j, TM1 [4l|] (which includes 
non- linear terms for the uj meson in order to soften the EOS at high densities), GM1 and 



GM3 



421, and FSU 



43j (with non linear u — p terms). For the second type we have 



considered the models TW 



44|, DD-ME1, DD-ME2 [45|L and DDH<5 [4Q]. The last one 



includes the 5-meson whose presence, as shown in Refs. 



471 ] and [481 ]. softens the symmetry 



energy at subsaturation densities and hardens it above saturation density. Finally, we have 
also considered the so-called quark-meson coupling model (QMC) [49]. In this model, nuclear 
matter is described as a system of non-overlapping MIT bags which interact through the 
exchange of scalar and vector mean fields. We have considered the parametrization used by 
Santos et al. in [50j, where saturation properties of asymmetric nuclear within QMC were 
compared with other relativistic models. Although having a quite high compressibility, the 
isovector channel properties such as the symmetry energy and its derivatives with respect 
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to the density are within the intervals set by experimental measurements. 



III. RESULTS 

We start the analysis of our results by showing in Tabled the bulk parameters character- 
izing the density dependence of the energy of symmetric matter and the symmetry energy 
around saturation density. Note that although all the models reasonably agree on their 
predictions for the energy per particle, symmetry energy, and density at saturation, they 
disagree in the remaining parameters, showing, in particular, significant discrepancies on 
the quantities Qo, L, K sym and Q sy m, and on the parameter K T = K sym — 6L — (Q /K )L 



that characterizes the isospin dependence of the incompressibility coefficient [56l. [57] ■ In the 
first three rows of the Table [J we show our BHF results with and without three-body forces 
respectively. As already mentioned the TBF contains two parameters that are fixed in order 
to reproduce the saturation point of symmetric matter. In the Table we present results for 



two sets of such parameters: the original set of Ref. 36J (labelled TBFa), and a second one 
(labelled TBFb), in which the parameter associated with the two pion attractive term has 
been reduced by 10%, and the one associated with the phenomenological repulsive term been 
increased by 20% in order to get a smaller saturation density. As a consequence, the binding 
energy and the incompressibility coefficient at the new saturation point decrease a little bit 
(see Table I). The slightly different results obtained with TBFa or with TBFb give an insight 
on the importance of the more phenomenological component of our approach. For most of 
the properties associated to the EoS, the differences are relatively small, which suggests that 
the results that we obtain are rather robust. The comparison of the different quantities is 
strongly influenced by the fact that they are calculated at the saturation density of each 
approach. Note that the effects of TBF are more important on the iso-scalar properties as 
Kq than on the properties associated to the density dependence of the symmetry energy. 
All the BHF results shown in the different figures contain the effects of TBF and have been 



obtained with the original set of parameters of Ref. 36] . 

Recent results from isospin diffusion (ID) predict values of L = 88 ± 25 MeV and K T = 
-500 ± 50 MeV Q, Q. The latter is in agreement with the value of K T = -550 ± 100 
MeV predicted by the independent measurement of the isotopic dependence of the giant 
monopole resonance (GMR) in Sn isotopes 0|. Similar values of L have also been 
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obtained by using experimental double ratios of proton and neutron spectra together with 



improved quantum molecular dynamics calculations 53]. While our Brueckner calculation 
leads to a value of L = 66.5 MeV, compatible with the data from isospin diffusion, our 
prediction of K T = —334.7 MeV is far from the lower bound K T = —450 MeV imposed 
by experiments. However, one has to be cautious when interpreting our parameter K T , 
defined in asymmetric nuclear matter, with the experimental data which, as pointed out by 
Blaizot and Grammaticos [54J], may be related to the isospin-dependent part of the surface 
properties of finite nuclei, especially the surface symmetry energy. Note that the effect of 
TBF on these two parameters is only of about 5 — 10%. Furthermore, only the Skyrme 
force SV and the relativistic models TM1, GM1 and QMC predict values of L and K T both 
compatible with experimental constrains, although their predictions for K are much larger 
than the value of Kq = 240 ± 20 MeV, supported nowadays by experimental data [55j . A 
close inspection to Tabled] shows, in fact, that none of the models considered predicts values 
of Kq, L and K T simultaneously consistent with t he p resent experimental data. In fact, 



as pointed out recently by several authors |56|, |57|, l58l . l59| |. it is difficult to determine the 



experimental value of K T accurately, since no single theoretical model is able to describe 
correctly the recent measurements of the isotopic dependence of GMR in Sn isotopes and, 
simultaneously, the GMR energies of a variety of nuclei. That suggests, as discussed by 



19 



52| 



Piekarewicz and Centelles in Ref. [56j, that the value of K T = —550 ± 100 MeV 
may suffer from the same ambiguities already encountered in earlier attempts to extract 
the incompressibility coefficient of infinite nuclear matter from finite-nuclei extrapolations. 
Concerning the third density derivative of the symmetry energy Q sym , for which there is not, 
at present, any experimental constraint, the microscopic prediction is large and negative 
which is in constrat with most of the effective models (except TM1) that give a large variety 
of positive values. 

In Fig. [3], we show the density dependence of S2G0) (left panel) for our BHF calculations 
and some representative Skyrme forces and relativistic models. In general, there is a good 
agreement between the microscopic BHF calculation and the Skyrme models considered in 
the whole density range explored. The relativistic models considered also agree with the BHF 
calculation at low densities, but deviations are found for TM1 and NL3 above saturation 
densities. A better insight of the density dependence of 5 2 (p) can be obtained by looking 



at the density dependence of the slope parameter L. The results for the same models are 
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plotted in the right panel of the figure. Note that even the models that apparently show a 
similar behaviour of the symmetry energy with density, as for instance the relativistic TW 
model and our BHF calculation, predict a different density dependence of L. In general the 
relativistic models predict a stiffer dependence of the symmetry energy, reflected in larger 
values of L than those produced in BHF. The Skyrme models shown in the figure, produce 
smaller values of L. These are selected models whose density behaviour of the symmetry 
energy has been tested [13]. 

We show in Fig. H]the correlations between L and K sym (left panel) and L and K T (right 
panel) , already considered in the literature for the case of effective forces [20|, |57| . Note that 
the BHF results for L are located inside the region delimited by the isospin diffusion data 
constraints, and that they adjust reasonably well to these correlations. There is no direct 
empirical information on the K sym parameter. However, as pointed out recently by Chen et 



al. in Ref. 



571 ]. whenever accurate experimental information becomes available for L, these 



correlations could be exploited to obtain theoretical estimates for K sym and K T . 

It has been shown by Brown and Typel 23|, 12J], and confirmed latter by other authors 



27 



43. 



51 



60 



61] , that the neutron skin thickness, 5R = J (r 2 ) — J(r%), calculated 



in mean field models with either non-relativistic or relativistic effective interactions is very 
sensitive to the density dependence of the nuclear symmetry energy, and, in particular, 
to the slope parameter L at the normal nuclear saturation density. Using our Brueckner 
approach and the several Skyrme forces and relativistic models considered we have made an 
estimation the neutron skin thickness of 208 Pb and 132 Sn and we have studied its correlation 
with the parameters L and K sym . However, a fully self-consistent finite nuclei calculation 
based on the BHF approach is too difficult to implement, therefore, following the suggestion 
of Steiner et al. in Ref. {3] we have estimated 5R to lowest order in the diffuseness corrections 
as SR ~ v^t, being t the thickness of semi-infinite asymmetric nuclear matter, 



t 



E s Jo 



Po(A0 1/2 



p^ 2 [E sym /S 2 (p) - l][E SNM (p) - E }^ 2 dp 



(15) 



Po(&)(l - P 2 c) 4vrr 2 f*>V°) p^[E SNM (p) - E ]^dp 

In the above expression E s is the surface energy taken from the semi-empirical mass formula 
equal to 17.23 MeV, r is obtained from the normalization condition (47rro/3)(0.16) = 1, and 
P c is the isospin asymmetry in the center of the nucleus. We have checked from Thomas- 
Fermi calculations that the value of (3 C is about 1/2 of the total isospin asymmetry of the 
nucleus /3 if Coulomb effects are not considered as in the present case. Therefore, we have 
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taken j3 c = (3/2. Although we can perform finite nuclei calculations with the Skyrme and 
relativistic effective models, we have also used Eq. f|T5|) with these models for consistency 
reasons with the BHF approach. We have checked that in the case of the Skyrme forces the 
accuracy of the results obtained by using Eq. ([T5]) with respect to a Hartree-Fock calculation 
is about 15 — 25%. 

In Fig. [5] we show the correlation between SR for 208 Pb (upper panels) and 132 Sn (lower 
panels) with the parameters L (left panels) and K sym (right panels). It can be seen, as 
has already been shown by other authors, that both the Skyrme forces and the relativistic 
models predict values of SR that exhibit a tight linear correlation with L. Note that the 
microscopic Brueckner calculation is in excellent agreement with this correlation. Since 
K sym is linearly correlated with L, it is, therefore, not surprising that 5R presents also an 
almost linear correlation with K sym , although it is less strong than its correlation with L, 
and shows a slight irregular behaviour. The linear increase of SR with L is not surprising 



since, as discussed in Refs 



27 



43 



51 



60 



61] , the thickness of the neutron skin 



in heavy nuclei is determined by the pressure difference between neutrons and protons, 
which is proportional to the parameter L, P(p ,j3) ~ Lp j3 2 /3 [20(. This can be seen for 
instance in Fig. [61 where we show the pressure of asymmetric nuclear matter evaluated at 
normal saturation density po for the isospin asymmetries of 208 Pb {(3 = 44/208) and 132 Sn 
(f3 = 32/132). Note that, since the BHF calculation predicts larger values for p than the 
other approaches (see Tabled]), its result appears in both cases, 208 Pb and 132 Sn, a bit out 
of the trends marked by the other models. 

Another sensitive quantity to the symmetry energy is the transition density pt from 
non-uniform to uniform /3-stable matter which may be estimated from the crossing of the 0- 



equilibrium equation of state with the thermodynamical spinodal instability line 



48 



64 



65. 



62. 



63, 



66j . As it has been shown in Ref . 65( the predictions for the transition density from 



the thermodynamical spinodal are ~ 15% larger than the value obtained from a Thomas- 
Fermi calculation of the pasta phase. Therefore, we may expect that our estimation of the 
transition density from the thermodynamical spinodal will define an upper bound to the 
true transition density |67j. We display in Fig. Ulpt as a function of the parameters L and 
Ksym for our BHF calculation together with the predictions of the several Skyrme forces and 
relativistic models. It is clear from the figure that p t is sensitive to the slope and curvature 
parameters L and K sym of the symmetry energy, decreasing almost linearly with increasing 
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L and K sym in agreement with recent results 66|, |68j. Using the experimental constraint 
on L from isospin diffusion, we estimate the transition density to be between 0.063 fm -3 
and 0.083 fm -3 . This range is in reasonable agreement with the the value of p t ~ 0.08 
fm -3 often used in the literature. Recently, Xu et al, 66] have obtained a different range 
for the transition density, namely, from 0.04 to 0.065 fm~ 3 using 51 Skyrme interactions. 
These authors argue that their results are smaller than 0.08 fm -3 because their approach is 
exact and no parabolic approximation is assumed for the isospin dependence of the nuclear 
force. However, we note that in the present work the parabolic assumption has only been 
considered in the BHF calculation whereas all the other results are exact and the obtained 
range of transition densities is in all cases the one indicated above. 

Finally, we show in Fig. [8] the transition density p t from non- uniform to /3-stable matter 
as a function of the neutron skin thickness in 208 Pb (left panel) and 132 Sn for our Brueckner 
calculation and the different Skyrme forces and relativistic models. The figure shows, as 
already pointed out by Horowitz and Piekarewicz 2^] that there is an inverse correlation 
between the neutron skin thickness and pt- In [27] a NLWM with non-linear u — p terms was 
used and the transition density was obtained with an RPA approach. We confirm the same 
trend for a larger set of nuclear models. Note that, again, our microscopic Brueckner results 
are in very good agreement with this correlation. As pointed out in Ref. 



271 ]. these results 

suggest that an accurate measurement of the neutron radius in heavy nuclei like 208 Pb or 
132 Sn is very important since it can provide considerable and valuable information on the 
thickness and other properties of neutron star's crust. 



IV. SUMMARY 



We have studied the density dependence of the symmetry energy within the microscopic 
Brueckner-Hartree-Fock approach using the realistic AV18 potential plus a three-body force 
of Urbana type. Our results have been compared with those obtained with several Skyrme 
forces and relativistic effective models. We have found a value of the slope parameter 
L = 66.5 MeV, compatible with recent experimental constraints from isospin diffusion, 
L = 88 ± 25 MeV. We have studied the correlation between the neutron skin thickness of 
neutron-rich isotopes and L and K sym . We have found that the BHF results are in very good 
agreement with the correlations already predicted by other authors using non-relativistic and 
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relativistic effective models. This agreement suggests that these correlations are not only due 
to the mean-field nature of these approaches. Microscopic calculations, as the one performed 
here, also provide a meaningful description of the isoscalar and isovector properties of the 
EoS and complement the already gathered knowledge on the bulk parameters of nuclear 
matter. We have also analyzed the correlations of L, K sym and the neutron skin thickness 
with the transition density pt from non-uniform to /3-stable matter in neutron stars. Using 
the experimental constraint on L from isospin diffusion, we have estimated the value of pt to 
be between 0.063 fm -3 and 0.083 fm -3 , a range in reasonable agreement with the the value 
of pt ~ 0.08 fm~ 3 often used in the literature. Finally, we have confirmed for a large set of 
nuclear models that there is an inverse correlation between the neutron skin thickness and 
the transition density p t , a trend pointed out first by Horowitz and Piekarewicz in Ref. j^ . 
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Model p E Kq Qo E sym L K sym Q sym K T Ref. 

BHF (with TBFa) 0.187 -15.23 195.5 -280.9 34.3 66.5 -31.3 -112.8 -334.7 

BHF (with TBFb) 0.176 -14.62 185.9 -224.9 33.6 66.9 -23.4 -162.8 -343.8 

BHF (without TBF) 0.240 -17.30 213.6 -225.1 35.8 63.1 -27.8 -159.8 -339.6 

SLy4 0.159 -15.97 229.8 -362.9 31.8 45.3 -119.8 520.8 -320.4 [8] 

SLylO 0.155 -15.90 229.7 -358.3 32.1 39.2 -142.4 590.9 -316.7 [9] 

SLy230a 0.160 -15.98 229.9 -364.2 31.8 43.9 -98.4 602.8 -292.7 [10] 

SkI4 0.162 -16.15 250.3 -335.7 29.6 59.9 -43.4 358.8 -322.5 [37] 

SkI5 0.156 -15.84 255.6 -301.7 36.4 128.9 159.8 11.2 -461.6 [37] 

SkI6 0.159 -15.88 248.2 -326.7 34.4 82.1 -0.9 332.3 -385.8 [38] 

SIII 0.145 -15.85 353.9 101.3 28.1 10.1 -392.3 130.2 -456.0 [39] 

SV 0.155 -16.04 305.3 -175.5 32.9 96.5 24.1 48.0 -499.4 [39] 

NL3 0.148 -16.24 271.6 203.1 37.4 118.5 100.9 181.2 -698.4 [40] 

TM1 0.145 -16.32 281.0 -285.2 36.8 110.8 33.6 -66.4 -518.7 [41] 

GM1 0.153 -16.32 299.2 -216.5 32.4 93.9 17.9 25.1 -477.5 [42] 

GM3 0.153 -16.32 239.7 -512.9 32.4 89.7 -6.5 55.8 -352.7 [42] 

FSU 0.148 -16.30 230.0 -523.4 32.6 60.5 -51.3 424.1 -276.6 [43J 

TW 0.153 -16.25 240.1 -540.1 32.7 55.3 -124.7 535.2 -332.1 [44] 

DDME1 0.152 -16.23 243.7 332.8 33.1 55.6 -100.8 703.8 -508.1 [45] 

DDME2 0.152 -16.14 250.8 478.1 32.3 51.4 -86.6 773.9 -493.8 [45J 

DDH<5 0.153 -16.25 240.2 -539.7 25.1 44.0 44.9 928.3 -120.2 [46] 

QMC 0.150 -15.70 291.0 -387.5 33.7 93.5 -10.0 28.0 -446.4 [49J 



TABLE I: Bulk parameters characterizing the density dependence of the energy of symmetric 
matter and the symmetry energy around saturation density for the BHF calculation with and 
without TBF and several Skyrme forces and relativistic models. All the quantities are in MeV, 
with the exception of po given in fm -3 . 
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FIG. 1: (Color online) Density dependence of the symmetry energy coefficients 52 and S4. Results 
for the BHF calculation, the Skyrme force SLy230a and the relativistic model NL3 are shown in 
the left, middle and right panels, respectively. 
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FIG. 2: (Color online) Isospin asymmetry dependece of the density (left panel), energy per particle 
(middle panel) and incompressibility coefficient (right panel) at the saturation point of asymmetric 
nuclear matter. Solid lines show the results of the exact BHF calculation whereas dashed lines 
indicate the results of the expansion of Eq. (fTUj) . Units of Eq((3) and Kq((3) are given in MeV 
whereas po{(3) is given in fm -3 . 
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FIG. 3: (Color online) Density dependence of the symmetry energy (left panel) and L (right panel) 
for the BHF calculation and some of the considered Skyrme forces and relativistic models. 
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FIG. 7: (Color online) Transition density from non-uniform to uniform /3-stable matter as a function 
of L (left panel) and K sym (right panel). The vertical dashed lines on the left panel denote the 
constraints on L from isospin diffusion (ID) 
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FIG. 8: (Color online) Transition density from non-uniform to uniform /3-stable matter versus the 
neutron skin thickness for 208 Pb (left panel) and 132 Sn (right panel). 
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